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Born-Infeld non-linear electrodynamics arises naturally as a field theory description of the dy- 
namics of strings and branes. Most analyses of this theory have been limited to studying it as 
a classical field theory. We quantize this theory on a Euclidean 4-dimensional space-time lattice 
and determine its properties using Monte-Carlo simulations. The electromagnetic field around a 
static point charge is measured using Liischer-Weisz methods to overcome the sign problem as- 
sociated with the introduction of this charge. The D field appears identical to that of Maxwell 
QED. However, the E field is enhanced by quantum fluctuations, while still showing the short dis- 
tance screening observed in the classical theory. In addition, whereas for the classical theory, the 
screening increases without bound as the non-linearity increases, the quantum theory approaches 
a limiting conformal fleld theory. 
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I. INTRODUCTION 



Theories of strings and branes have shown much promise as theories unifying gravity with 
strong, weak and electromagnetic interactions. One hmitation has been that most of the 
work on these theories has been performed on the classical theories. The most promising 
way to produce quantum theories of strings and branes has been to study field theories 
which describe their dynamics. The field theories which arise most naturally in this context 
are Born-Infeld non-linear electrodynamics and its extensions. Early papers which use 
Born-Infeld theories to describe the dynamics of strings and branes include 

mm. 

However, since much of the interesting physics arising from the non-linearities in Born-Infeld 
models is inaccessible to perturbation theory, most of the work on these theories has also 



been limited to the classical domain 



m 



or to the small quantum fluctuations 9] around 
these classical solutions. Beyond this, there have only been a few exploratory investigations 
of the general problem of quantizing this theory, for example , which examines quantum 
solutions representing propagation in a fixed direction. 

Born-Infeld electrodynamics in + 1 dimensions is described by the non-linear action 



(1) 



To make the connection with strings and branes, one typically chooses n = 9 and dimen- 
sionally reduces the theory from 9 + 1 down to p + 1 dimensions. This describes a p-brane, 
and the extra 9 — p components of are identified with the transverse degrees of freedom 
of the brane. The special case p = 1 describes a string. Here the non-linearity parameter h 
is related to the string tension through 1/6 = 2Tia' . 

In this work we study the original Born-Infeld theory in 3 + 1 dimensions [n = 3) with no 
dimensional reduction. In brane language we study a 3-brane with no transverse dimensions. 
This theory was originally proposed as a modification of QED in which the electric field of 
a static point particle was screened at short distances, rendering the self-energy of a point 
particle finite [l[Q]. 

We note that, in Euclidean space, the action is positive so that quantization by the 
standard functional integral techniques is well defined. We then define this theory on a 
discrete space-time lattice, preserving gauge invariance. The positivity of the action allows 
the use of Metropolis Monte-Carlo JJj simulations to extract the properties of the quantum 
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field theory. In particular, we measure the electromagnetic fields produced by a static point 
charge by including the eff'ects of a Wilson Line (Polyakov Loop) in the action. This destroys 
the positivity, but we are able to use the methods of Liischer and Weisz (after Parisi, 
Petronzio and Rapuano 13]) to produce the exponential statistics required to overcome this 
lack of positivity. 

We have performed simulations of the lattice Born-Infeld theory on 8^ and 12^ lattices. 
We present evidence that the D field emanating from a point charge is identical with that 
of the Maxwell theory, for all values of b, up to lattice artifacts. The E field shows short- 
distance screening as for the classical theory, but whereas classically the screening length 
increases without limit as 6 ^ 0, the quantum field theory approaches a conformal field 
theory in this limit. In addition, the E field is enhanced over its classical value by quantum 
fiuctuations. Preliminary results for the smaller lattice were presented at Lattice 2005, 
Dublin Q- 

In section 2 we summarise the salient features of the classical theory from the extensive 
literature on Born-Infeld electrodynamics. Section 3 discusses the lattice formulation in 
Euclidean space and simulation methods. We present the results of our simulations in 
section 4. These results are discussed and conclusions drawn in section 5. 



II. CLASSICAL BORN-INFELD ELECTRODYNAMICS 



In this section we summarise those results from the literature on the classical Born- 
Infeld theory in 3 + 1 dimensional Minkowski space which are relevant to our investigations. 
These results are condensed from the following references 
determinant in equation ^ leads to the Lagrangian 
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171, Q- Evaluating the 
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1-Jl- b-^{E^ - B2) - b-^{E • B) 



(2) 



Note that the requirement that the argument of the square root be positive, restricts the 
magnitude of E which in turn gives rise to short-distance screening and leads to a finite 
self energy for a point charge. This was the original reason Born and Infeld proposed this 
non-linear extension of electrodynamics. When dimensionally reduced Born-Infeld theories 
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are used to describe strings or branes, this same mechanism restricts the energy density of 

the string/brane to be finite. 

To proceed, it is customary to define D and H fields by 

_ dC _ E + b-\E • B)B 

~ dE~ li _ 5-2(^2 _ B^) - h-HE ■ BV 



H = —= B-h-\E-B)E 

dB . Il _ 5-2(^2 _ _B2) _ ^-4(^ . ' 



In terms of E, D, B and H, the field equations are the standard Maxwell equations 

V -D = p 

V -B = 

^ ^ dB 

V X E+— = 

ot 

dD 

V.H--^j, (4) 

where, as usual, the interaction with the electromagnetic current is introduced by adding 
a term j^A^j^ to the Lagrangian. All the non-linearity is hidden in the defining relations of 
equations El 

For a static point charge p = e6^{r), V • D = p has the spherically symmetric solution 

D = -^r (5) 
as for the Maxwell case. B = H = and equations El then yield 

where 



defines the screening length. Due to this short distance screening, the electric field at the 
origin has magnitude b rather than being infinite, and the energy in the field of the point 
charge is finite. 

III. EUCLIDEAN BORN-INFELD ELECTRODYNAMICS ON A LATTICE 

In Euclidean space, the Born-Infeld action is 



S = d^x 
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(8) 



We note that this action is positive, so this theory can be simulated using importance 
samphng methods. The expressions for D and H are identical to those of equations 01 
except that all signs are positive. 

To enable simulations, we must first transcribe this field theory to a discrete space-time 
lattice. We choose a hypercubic lattice with lattice spacing a. For convenience we define 
(3 = b'^a^. From now on we will work in lattice units where a is set to 1, except where 
we discuss the effects of varying a. Because power counting suggests this theory is non- 
renormalizable, it is generally considered as an effective theory, requiring a cutoff. Hence we 
are not required to take the a — > limit. We shall have more to say about this later. 

We choose the non-compact formulation of QED on the lattice, since this is closest to 
continuum QED (In the Maxwell case it defines a solvable, free field theory, as in the 
continuum). This is done by defining the lattice F^^^, by 



which is gauge invariant as is the continuum F^^,. This indicates the first subtlety; F^,y is 
not associated with a single site x or link of the original lattice. Therefore to preserve the 
symmetries of the cubic lattice we average our action over the 16 choices of 6 plaquettes 
associated with a lattice site. These 16 choices are defined by choosing 1 unit vector in 
each of the 4 directions emanating from the site, in all possible ways. The set of plaquettes 
associated with a given choice of unit vectors is uniquely defined by requiring that 2 edges 
of each plaquette belong to the chosen set of unit vectors. Simulations are performed using 
Metropolis Monte-Carlo updates of the gauge fields (A^) on each of the links of the lattice. 
We periodically perform subtractions of constant fields from each A^, since such constant 
fields have no physics, to prevent the lattice average of A^ becoming too large. Similarly, 
we periodically gauge-fix these fields to Landau gauge, to remove large gauge fiuctuations. 

To compare with the classical Born-Infeld results we need to measure the E and D fields 
for a point charge. Since we use periodic boundary conditions on the gauge fields, the 
total charge on the lattice must be zero. We choose the 'Jellium' solution from condensed 
matter physics to circumvent this difficulty. Here, if the magnitude of the point charge is 
e, we distribute a charge — e uniformly over the lattice. Hence each site has charge — e/V 



F^u{x + l/t + |z>) = Au{x + ft) ~ Au{x) - A^{x + z>) + A^{x), 



(9) 
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(V = NxNyNz is the spatial volume of the lattice), except for that site containing the point 
charge, whose charge is now reduced to e — e/V. On the lattice this charge is introduced by 
including a Wilson Line (Polyakov loop) W{x) in the functional integral. 



W{x) = exp < ie 



''^ y 



(10) 



We calculate the {E) and (D) due to this charge as 

{E{y,t)W{x)) 



i{E)piy-x) 



{Wix)) 



where the i converts our electric fields from Euclidean space to Minkowski space, and we 
use the subscript p to denote the expectation value in the presence of charge density p. 
The Wilson Line for fixed x on any given configuration has magnitude one, so 

{\W{x)\) = l, (12) 

while it is expected that 

(Wix)) = e-"'"'''^'. (13) 

This indicates that the phase of W causes important cancellations between the Wilson lines 
for different configurations in our ensemble. Hence we have a sign problem. Averaging over 
the site x helps improve our statistics. Even if we measured the Wilson line averaged over 
all sites, every sweep, we will run into trouble if 



{W{x)) ^ l/y (number of sweeps x V, (14) 

in the most optimistic scenario where the phases of the Wilson lines on a given configuration 
are uncorrelated as are those on consecutive sweeps. For any of our 12^ runs, the number 
of sweeps is 500,000 and V = 1728 and the right hand side of equation IT^ is ~ 3.4 x 10~^. 
As we shall see later, we really need to do better. 

Liischer and Weisz jl^ (following Parisi, Petronzio and Rapuano 'l^ ) have pointed out 
that one can effectively enhance ones statistics by dividing the lattice into some number (call 
it n) of timeslices. Notice that, if one fixes the link fields on the boundaries of these timeslices, 
the Monte Carlo updates of those fields on the interiors of each slice are independent of the 
updates of the interiors of all other slices. Thus if one performs m updates of the interior 



of each time-slice, it is equivalent to performing m" (restricted) updates of the lattice. In 
addition this process can be performed recursively over progressively coarser slicings of the 
lattice. In practice we have used thickness 1 and 2 time slices recursively. 

Let us now examine the limiting cases, first where b oo, and second where 6-^0. As 
b ^ oo, the Euclidean Lagrangian approaches its Maxwell (free field) form 

Cj, = ^{E' + B'). (15) 

In the limit b ^ the Euclidean Lagrangian approaches the limiting form 

Ce = \E-B\, (16) 

which describes a 4-dimensional conformal field theory. [Note that the Minkowski La- 
grangian (equation P) does not exist in this limit. However, the Hamiltonian formulation is 
valid in this limit, and leads to the Hamiltonian 

n=\DxB\, (17) 

which also describes a conformal field theory 
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IV. SIMULATION DETAILS AND RESULTS 

We have simulated lattice Born-Infeld QED on 8^ and 12^ lattices. On both 8^ and 
12^ lattices we performed runs of 500,000 10-hit Metropolis sweeps of the lattice at each of 
f3 = 0.0001,0.01, 1, 100. On the 8'^ lattice we ran an additional 100,000 sweeps at each of 
f3 = 0.1,0.2,0.5,2,5. Measurements were made every 100 sweeps. On the 8^ lattices the 
Wilson lines and the on-axis fields they produced were measured in the Liischer-Weisz scheme 
where we used 10 10-hit Metropolis updates of the interiors of each thickness-1 timeslice for 
each of 10 10-hit Metropolis updates of the interiors of each thickness-2 timeslice. For the 12'^ 
lattice we increased the number of updates at each level to 20. This means that we effectively 
used 10^^ configurations every 100 sweeps for our 8^ lattice simulations and 2.62144 x 10^^ 
configurations every 100 sweeps for our 12^ runs, for our measurements of the Wilson lines 
and electric fields. Admittedly these configurations generated in the Liischer-Weisz scheme 
are not all independent. 

In figure[T]we present the average action density {S) /V, where V is the space-time volume, 
as a function of (3 = b'^a^ for our runs on both 8^ and 12^ lattices. We notice that there is 
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FIG. 1: Euclidean action density {S) jV as a function of /3 

good agreement between the 2 lattice sizes, indicating that that the finite volume effects are 
small. As /3 ^ oo, where the theory becomes the linear Maxwell theory, the equipartition 
theorem predicts that the action density is |. We note that the action density is within 
a percent of this value by /? = 100, indicating that for /5 ^ 100 the theory is close to the 
Maxwell theory. In the /5 — > limit the scaling properties of the action again lead to an 
action density of |. We see that our = 0.0001 action density is very close to this value, 
while the f3 = 0.01 value is only a little over 2% higher. This suggests that lattice Born-Infeld 
QED is close to the limit described by equation Uni for j3 ^ 0.01. [Note that, because the zero 
modes do not contribute to the action, the limiting values (|) have corrections (9(1/V)]. In 
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between the 2 limits the action rises to a maximum of magnitude < 2 for /3 ~ 1. 

Figure |21 shows the Wilson Lines as functions of electric charge e for 16 values of this 
charge ranging from 0.0078125 to 4.0 and for the 4 beta values /? = 0.0001, 0.01, 1.0, 100.0. 
At each (3, (W) starts from 1 at e = 0, and falls faster than exponentially with charge e as e 
is increased. The falloff becomes faster as (3 (or b) is decreased, and the field theory becomes 
more non-linear. We note, however, that the values of the Wilson lines for (3 = 0.0001 
and (3 = 0.01 are almost identical. This is further evidence that the theory is close to the 
f3 = theory for f3 ^ 0.01. This also gives us the first evidence for differences between the 
electrostatics of the quantum and classical theories. Classical electrostatics has B and H 
identically zero, the contribution of equation (TBI vanishes and the argument the square root 
of equation n reduces to 1 — b~'^E^ for all b. In the quantum theory it is not even consistent 
to set B = 0, and the result is that the action still reduces to the independent form 
of equation Uni as 6 — » 0. We note that the Liischer-Weisz prescription yields an excellent 
signal for the Wilson line with small error bars over almost 11 orders of magnitude. 

We measure the E field for on axis separations (for definiteness consider this to be in 
the z direction) from the electric charge. Here we need only measure Ez- For this we 
need to measure plaquettes in the {z, t) plane containing the Wilson Line. We call the z 
separation of the centre of such a plaquette from the Wilson line Z. Clearly Z takes the 
values |, |, |, Nz — |. It is also clear that 

{Ez{Z)) = -{Ez{Nz-Z)) 

{Dz{Z)) = -{Dz{Nz-Z)). (18) 

The lattice field is obtained by calculating where E^ is a {z,t) plaquette and S is 
the lattice action. The {Ez)p and (-Dz)p are obtained using equationHH Again, the Liischer- 
Weisz method yields a good signal. In what follows, we will drop the subscripts p and z for 
convenience. 

We will start by considering D. Since the equation of motion 'V-D = p holds for all /5, this 
means that, to the extent that we have rotational invariance, (D) should be independent 
of /3, and equal to the classical value. Of course, since we are on the lattice, rotational 
invariance is at best approximate, and we expect departures from this expectation. 

Figure El shows the ratios of the D fields produced by a point charge e on both 8'* and 
12"^ lattices for /3 = 0.0001, where we expect the action to closely approximate the small {3 
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FIG. 2: Wilson Lines as functions of charge e, for a range of (3 = a%'^ values: a) 8^ lattice: b) 12"^ 
lattice. 
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FIG. 3: Ratios of D to their Maxwell (free field) values for /3 = 0.0001 as functions of electric 
charge e, for accessible values of the on-axis separation Z: a) for an 8^ lattice, b) for a 12^ lattice. 
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limit of equation to their free field values. For the minimum separation of source and 
field [Z = ^) D is consistent with its free-field (Maxwell) value, evaluated on the same-size 
lattice. This is because, at this separation, the discrete Gauss' equation plus cubic symmetry 
is sufficient to uniquely determine D; in fact D = |(1 — 1/V)e. At separation Z = D 
has fallen to 86-87% of its free-field value. Gauss' law and cubic symmetry are no longer 
sufficient to determine the D field uniquely, and we do not have the simple relation D = E 
and hence V x + dB /dt = of the Maxwell theory to help us. It is thus reassuring that 
D is as close to its free-field value as it is. As the separation increases, the ratio of D to its 
free- field value appears to be approaching 1. This is an indication that rotational symmetry 
is restored at length scales much larger than the lattice spacing. The fact that the D for 
this highly non-linear theory appears identical to that of the Maxwell theory, apart from 
lattice artifacts, is non-trivial. For the Maxwell theory, D = E, while in the (3—^0 limit 
D = e{E • B)B, where e is the familiar sign function. 

Figure m shows the e dependence of D, for Z < Nz/2, for P = 1.0. This (3 is close to that 
which maximizes (5*) /V. Thus, for this value of (3, Born-Infeld QED should be far from 
regular (Maxwell) QED, and also far from the non-linear limit of equation ^1 Again we 
see that at minimum separation, D is at its free field value. The ratio of D to its free-field 
value falls to a minimum for 1.5 ^ Z ^ 2.5, after which it begins to increase again. For the 
8^ lattice, the finite lattice extent tends to suppress this increase while the larger 12^ lattice 
shows a clear increase. Hence we see evidence that rotational symmetry is being restored at 
large distances, just as for small /3. Again we believe that the difference between the /? = 1.0 
D fields and their free-field counterparts is a lattice artifact. 

In figure El we graph the ratios of D to its free field value for f3 = 100.0, where we expect 
Born-Infeld QED to closely approximate standard (Maxwell) QED. We note again that at 
minimum separation D is consistent with its free field value. At larger separations there are 
small but significant departures from the Maxwell theory. However, since these departures 
are always less than (and often considerably less than) one percent, this is further evidence 
that (3 = 100.0 Born-Infeld QED is close to standard QED. Because these departures com- 
parable in size with our error-bars, it is difficult to make any more quantitative observations, 
although it does appear that the 8'^ and 12^ values are consistent for the ranges of Z which 
they share. 

We now turn our attention to the E field. Here, unlike the case of the D field. Maxwell's 
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FIG. 4: Ratios of D to their Maxwell (free field) values for /3 = 1.0 as functions of electric charge 
e, for accessible values of the on-axis separation Z: a) for an 8^ lattice, b) for a 12^ lattice. 
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equations are of little help, because the relations of equations El are highly non-linear. What 
we are interested in knowing is how similar these E fields are to their classical counterparts 
and, in particular, whether they show the short-distance screening of classical Born-Infeld 
electrodynamics. We shall use our experience with the D fields as a guide to how much the 
results that we observe are affected by lattice artifacts. 

First we look at the case where (3 = 0.0001 which we expect to be close to the non-linear 
limit of equation Figure IHl shows the ratios of E to its free-field value for this (3 as 
functions of e for all accessible values of Z. The similar graph (not shown) for f] = 0.01 
is almost identical, which justifies our claim that we are close to the non-linear limit. We 
first note that the 8'^ and 12^ graphs are very similar for those Z values common to each, 
indicating that we should be able to extract information which is not dominated by finite 
volume effects from lattices of this size. Secondly we notice that for all accessible Zs and 
over most of the range of charges considered (0 < e < 3), the electric fields are larger than 
their free field counterparts by a factor of 2-3. This contrasts with the classical case where 
E < D. This enhancement of E is due to quantum fluctuations which are clearly quite large. 
At Z = |, it is clear that E is screened, and that this screening increases with increasing 
e. The Z = ^ field also shows such screening, but the falloff with increasing e is much less 
pronounced than at Z = |, while there is no clear signal for screening for larger Z. This 
behaviour is qualitatively similar to what is observed classically. However, classically its 
screening is determined by the ratio e/b (or tq), whereas in the quantum theory, for small 
enough b (and hence /3) the screening becomes b independent. Since the limit (3^0 can be 
achieved by taking a at fixed b, rather than b —>■ 0, the non-linear limit of equation ITHl 
could be expected to describe the continuum limit. The question arises as to whether this 
limiting field theory is trivial or not. If it is trivial, the E field of a point charge (a propagator 
of this field theory) would be proportional to the free field propagator. Of course, this would 
only be true of this propagator at finite distances in physical units, which as a — means 
large lattice distances. Thus the fact that this ratio on the 12^ lattice varies by only of order 
5% for Z > 0.5, at least for small charges, could be an indication that this limit is trivial. 

Figure [7| shows E/Efree as functions of e for all accessible separations Z, for j3 = 1.0. 
Again, we note the similarity between the results from the 8^ and 12^ lattices up until Z 
approaches A^z/2 where the effect of the finite volume becomes apparent, which gives us 
confidence that the finite volume effects are under control. Despite the fact that we are now 
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FIG. 7: Ratios of E to their Maxwell (free field) values for (3 = 1.0 as functions of electric charge 
e, for accessible values of the on-axis separation Z: a) for an 8^ lattice, b) for a 12^ lattice. 
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far from both the small and large (3 limits, these graphs are qualitatively similar to those 
for (3 = 0.0001, except for the scale on the vertical axis. The enhancement of this electric 
field over its Maxwell value due to quantum fluctuations is smaller than that for (3 = 0.0001. 



over the same range of charges e. 

In the last of this series of figures showing the e and Z dependence of the ratio of E to 
its free field value, figure |Hl we show the results for (3 = 100 where we expect Born-Infeld 
QED to be close to its Maxwell (standard QED) limit. What is clear from these graphs is 
that, since the E never differs from its free field value by more than 2% over the range of 
e considered, Born-Infeld QED is close to its free-field limit for (3 = 100. We do, however, 
observe a small increase in E over its free-field value due to quantum fluctuations, and see 
clear evidence for weak screening for Z = ^. This is good evidence that expanding the 
square root in the action in powers of 1//5, and using a perturbative analysis, would be valid 
in this domain. 

In an effort to get a more quantitative understanding of the screening of the E field that 
we have observed, we look again at the screening in classical Born-Infeld electrodynamics. 
Here we can expose the screening length (rg) by looking at 



So, for the classical case, plotting the left hand side of this equation for fixed r against 
(or I/tq) would give a straight line. Figure IHl shows this quantity from our simulations at 



we have drawn to guide the eye indicate that the dependence is almost linear. Perhaps 
the observed departures are due to lattice artifacts. What this behaviour does suggest is 
that the screening increases without bound as is increased, as is true classically. 

Finally, for completeness, we show the Wilson Line correlations for each of the 16 charges 
for P = 100 and (3 = 0.0001 on the 12'^ lattice in figure ^1 These correlation functions 
show the interactions between charge +e and charge — e separated by distance Z. What 
these indicate is that the effect of screening on such interactions is to reduce the correlation 
functions as (3 is decreased (which as observed above, increases screening). One expects that 
the Wilson Line correlation function will behave as 



The screening which we observe at Z = ^ (and Z = |) is also less by roughly a factor of 3 




(19) 



j3 = 0.0001. While straight line fits to the Z = and Z = ^ are impossible, the straight lines 



{W{x)W^{x + Z)) = C{Z) exp[-V(Z)T]. 



(20) 
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Although the electric field shows screening, E for /3 — 0.0001 exceeds that for the free field 
and hence that for /3 = 100.0 for the same charge, over the range of charges considered. 
Hence, the potential energy V{Z) at P — 0.0001 exceeds that for /3 = 100.0. It is for this 
reason that this correlation function falls more rapidly for the smaller /?. 

V. DISCUSSION AND CONCLUSIONS 

Euclidean Born-Infeld electrodynamics has been quantized on a discrete space-time lat- 
tice. The Euclidean action is positive, allowing us to simulate this theory by importance 
sampling techniques. We have simulated the gauge-invariant non-compact lattice implemen- 
tation on 8^ and 12^ lattices using standard Metropolis Monte-Carlo methods for several 
values of the non-linearity parameter b [or the dimensionless j3 = h'^a^ (a is the lattice spac- 
ing)]. We measure the electric fields produced by a static point charge e, as functions of that 
charge. On the lattice this charge is introduced by the inclusion of a Wilson line (Polyakov 
loop). The phase of this Wilson line introduces a sign problem. However, this sign problem 
is overcome using the method of Liischer and Weisz, which uses partial factorization of the 
functional integral to produce the exponential statistics required to circumvent such sign 
problems. With this, wc were able to obtain excellent signals for the Wilson line, E and D 
fields, over an appreciable range of e values. 

We performed simulations with (3 values ranging from those large enough to approximate 
the Maxwell (free-field) theory to the those where Born-Infeld QED is well approximated 
by a /3(6)-independent conformal field theory with Euchdean Lagrangian Ce — \E • B\ or 
(in Minkowski space) by the Hamiltonian = |Z> x B|. For all (5 values, the electric dis- 
placement D appears to differ from its value in the Maxwell theory only by lattice artifacts. 
The electric field E shows considerable enhancement over its Maxwell value by as much as 
a factor of ~ 3 due to quantum fiuctuations. Away from the Maxwell limit, E exhibits 
short-distance screening with a screening length whose 4th power is approximately linear in 
e^, as in the classical theory. Also as in the classical theory, screening is enhanced as /3(6) 
is increased. However, unlike in the classical case where the screening at fixed e increases 
without limit as h is increased, in the quantum theory this screening approaches a limit 
described by the conformal field theory mentioned at the beginning of this paragraph. Once 
this limit is achieved, the only way to increase the screening length is to increase |e|. 
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Since our chief reason for studying Born-Infeld QED is because of its connection to strings 
and branes, we only need to consider it as an effective field theory with a momentum cutoff, 
or in our language, a finite lattice spacing. However, it is interesting to consider whether one 
can remove the momentum cutoff (take a ^ 0) in such a way as to define a rcnormalizable 
field theory. If we fix b and take a to zero, we are taking the (3^0 limit, where the field 
theory is described by the Euclidean Lagrangian Ce — \E • B\. Hence this theory describes 
the infinite-momentum-cutoff limit of Born-Infeld QED. What we do not know, is whether 
this quantum field theory is non-trivial. Our simulations hint that the propagators of this 
theory might show free-field scaling. However, the evidence is far from convincing, especially 
if we consider the possibility that deviations from free-field scaling could well be logarithmic. 

It would be useful to perform simulations of Born-Infeld QED on larger lattices, possibly 
restricting these simulations to those using the action based on the Lagrangian Ce — \E'B\. 
Of more interest would be simulations using those actions obtained from higher dimensional 
Born-Infeld theories by dimensional reduction, which describe strings or branes with trans- 
verse degrees of freedom. In these actions, the square root contains an additional term 
{d^Xa^, where the scalar fields Xa are the n — p additional components of the A field cor- 
responding to the reduced dimensions and are interpreted as the transverse displacements 
of the p-brane. We plan such simulations in the near future. 
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